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Abstract. We generalize the results of Song-Zelditch on geodesies in spaces of 
Kahler metrics on toric varieties to harmonic maps of any compact Riemannian 
manifold with boundary into the space of Kahler metrics on a toric variety. We 
show that the harmonic map equation can always be solved and that such maps 
may be approximated in the C 2 topology by harmonic maps into the spaces 
of Bergman metrics. In particular, WZW maps, or equivalently solutions of a 
homogeneous Monge- Ampere equation on the product of the manifold with a 
Riemann surface with S 1 boundary admit such approximations. We also show 
that the Eells-Sampson flow on the space of Kahler potentials is transformed to 
the usual heat flow on the space of symplectic potentials under the Legendre 
transform. 



1. Introduction 

Our main purpose in this article is to prove that the Dirichlet problem for a 
harmonic map tp : N —* 7i(T) of any compact Riemannian manifold with boundary 
N into the infinite-dimensional space 7i(T) of toric Kahler metrics on a smooth 
projective toric variety (M, oj) admits a smooth solution that may be approximated 
in C 2 (N x M) by a special squence of harmonic maps ipk : N — > Bk{T) C 7i(T) 
into the finite-dimensional subspaces of Bergman (or Fubini-Study) metrics induced 
from projective embeddings. As a special case, we show that the WZW (Wess- 
Zumino-Witten) equation, or equivalently the homogeneous complex Monge- Ampere 
equation on the product of the manifold with a Riemann surface with S 1 boundary, 
admits such a solution as well as such approximations. This generalizes previous 
work of Song-Zelditch in the case of geodesies, i.e., where N = [0, 1]. 

Before stating our results, we briefly recall the background to our problem. Let 
[M , u) be a compact closed Kahler manifold of dimension m with integral Kahler 
form and let (L,ho) — * M be an ample Hermitian holomorphic line bundle with 
Uh = uj satisfying [u] = ci(L), where = —^^ddlogho is the curvature (1,1)- 
form of ho. Any other hermitian metric on L may be expressed as hm = e~ tp hQ, with 
ip a smooth function on M. Following Mabuchi, Semmes and Donaldson [M, S2| ID1| 
one may regard the space 
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:= W e C°°(M) :u v = uj + -^dBif > } (1) 

Z7T 

of potentials of Kahler metrics in a fixed cohomology class as an infinite-dimensional 
symmetric space dual to the group of Hamiltonian diffeomorphisms Ham(M, uj). We 
will henceforth usually identify H u with the space of Hermitian metrics on L of 
positive curvature 



H:={h : h = h e- v ,tpe H^}. 
The symmetric space Riemannian metric g^2 is defined by 

9ip(C,ri) v :=^ f Cr?<, (,V ^ TpHui — C°°(M). (2) 

A basic idea, considered by Yau, Tian and Donaldson, is to approximate transcen- 
dental objects defined on TL by algebraic objects defined on the finite dimensional 
symmetric spaces Bk of Bergman (or Fubini-Study) metrics on L. To define them, 
we use the following notation: H°(M,L k ) is the space of holomorphic sections of 
the fc-th power L k -> M, d k + 1 = dimF°(M, L k ) and BH°(M, L k ) is the manifold 
of all bases s = {sq, . . . , Sd k } of H°(M, L k ). A basis s determines a Bergman metric 

h s := (4%s) 1/fe = -TT^r, (3) 

(E^I-iWIJ*) 

as the pullback of the Fubini-Study metric hps on the hyperplane bundle O(l) — ► ¥ dk 
under the Kodaira embedding 

Ls-.M^V 1 *, z^[s (z),...,s dk (z)}. (4) 

The space of all Bergman metrics defined by BH°(M, L k ) is denoted by 

B k = {hs, s€BH°(M,L k )}. (5) 

We may identify the space Bk with the symmetric space GL{dk + l,C)/U(dj i + 1) 
since GL{dk + 1,C) acts transitively on the set of bases, while tJ/iFS is unchanged 
if we replace the basis s by a unitary change of basis. 

A natural question is: to what extent can the geometry of Ti be approximated 
by that of the spaces Bk of 'algebro-geometric' metrics? At the most basic level 
of individual points, it follows by the Tian Asymptotic Isometry Theorem [Tj . and 
its subsequent refinements (Cj [Zj, that a metric h G TL can be approximated in a 
canonical way by a sequence of Bergman metrics hk ■ The geometry of a Riemannian 
manifold is reflected to a large extent by its geodesies and more generally by the 
specification of the harmonic maps into it, involving the analysis of certain nonlinear 
elliptic PDEs. In this article we describe how solutions to these PDEs on TL can be 
approximated in an algebro-geometric manner by a sequence of solutions to PDEs 
on Bj., in the setting of a toric variety. 
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We will need the canonical sequence mentioned above to state our results, so we 
recall how it is constructed. First, observe that Bk is isomorphic to the symmet- 
ric space Ik of Hermitian inner products on H°(M,L k ), the correspondence being 
that a basis is identified with an inner product for which the basis is Hermitian 
orthonormal. Define the maps 

Hilb fc : TL - l k (6) 

by the rule that a Hermitian metric h G TL induces the metrics h k on L k and the 
inner products on H°(M,L k ), 

where V = § m lo™. An inner product I = (•, •) on H°(M,L k ) determines an I- 
orthonormal basis s = Sj of H°(M, L k ), an associated Kodaira embedding as 
well as a Bergman metric given by 

FSfc(I) := ^. (8) 

Tian's asymptotic isometry theorem then states that FS& o Hilb^/i) — > h in the 
C°° (M) topology with complete asymptotic expansions [T l [C | [Z]. 

The question was then raised \AT\ \D2\ IPS1| whether geodesies of TL could be 
well-approximated by the one-parameter subgroup geodesies of Bk- Geodesies of TL 
are given by solutions of 

tp - ^\Vtp\ 2 = 0, ip(0, •) = ipo, tp(l, •) = ipi, ipo^xGH, (9) 

where cp is considered as a map from [0, 1] to Ti, or equivalently as a function on 
[0, 1] x M. Phong-Sturm studied this problem in depth and proved that a sequence 
of geodesies in Bk converges weakly (almost everywhere) to a prescribed geodesic of 
H |PSlj . Song-Zelditch proved that the same sequence converges in C 2 ([0, 1] x M) 
when the manifold is toric and one restricts to the torus-invariant metrics [SoZ2j, 
and Berndtsson used a different argument to prove that geodesies in TL can be C°- 
approximated by geodesies in spaces of Bergman metrics induced by embeddings by 
sections of L k ®KM, where Km is the canonical bundle of M [B]. In addition, Phong- 
Sturm and Song-Zelditch proved approximation results for geodesic rays constructed 
from test configurations [PS2, SoZ3]. 

A harmonic map between two Riemannian manifolds (N, f) and (iV, /) is a critical 
point of the energy functional 

E(a) = J \da\ 2 fm *f dv N,f, 

on the space of smooth maps a from N to N |ES] . The problem we study in this 
article is whether higher dimensional harmonic maps of general compact Riemannian 
manifolds N with boundary dN into TL admit similar kinds of 'algebro-geometric' 
approximations. For maps into the space of toric Kahler metrics on toric manifolds, 
we obtain an affirmative solution at the same level of precision as in the case of 
geodesies studied by Song-Zelditch. 
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To describe our results, let us define our setting more precisely. We recall that a 
toric variety M of complex dimension m carries the holomorphic action of a complex 
torus (C*) m with an open dense orbit. We let T = (S 1 ) 171 be the associated real 
torus. Objects associated to M are called toric if they are invariant with respect to 
T. We let uj denote a toric integral Kahler form on M, and let L be an ample line 
bundle with [u] = c\(L). We then define the space of toric Hermitian metrics on L, 



H(T) ={h£H:t*h = h, Vt £ T}. (10) 

This is a flat submanifold of TC. As before we will frequently identify an element of 
TC(T) with the corresponding element oi7i u . Moreover, we will often identify an ele- 
ment with the local Kahler potential defined on the open orbit of the complex torus; 
see Section El We also denote by Bk (T) C Bk the subspace of Bergman metrics de- 
fined by T invariant inner products, or equivalently by T-equivariant embeddings. 
Any such embedding is induced by the basis of toric monomials {x«}aefcPnz m 

of 

H°(M, L k ), where P denotes the moment polytope associated to the action, and 
kP denotes its dilation by a factor of k. Finally, let (N, f) be a compact oriented 
Riemannian manifold with smooth boundary, let G(y, q) denote the positive Dirich- 
let Green kernel for the Laplacian A^r := &N,f ( see Section and let dVgNj denote 
the induced measure on dN from the restriction of the Riemannian volume form 
dVpfj from N to dN. The main result of this article is: 

Theorem 1.1. Let (M,L,co) be a polarized toric Kahler manifold, and let (N,f) 
be a compact oriented smooth Riemannian manifold with smooth boundary dN. Let 
if) : dN — > 7i(T) denote a fixed smooth map. There exists a harmonic map <p : 
N — > TL{T) with ip\dN = and harmonic maps (pj~ : N — * Bk(T) with ^k\dN = 
FSfc o Hilbfe(^), given on the open orbit by 

<Pk{y,z) = T-log V |XaO)|j!fc exp C / d u G(y,q)log\\xa\\lk dV dNJ (q)), 
a&kPHZ™ JdN 
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and one has 



in the C 2 (N x M) topology. 



lim (p k = ip, 

k— >oo 



A motivating special case is the unit disc N = D := {z £ C : \z\ < 1}. It 
has been the subject of intensive studies (e.g., \Ch\ ICT[ ID3] ). Then the map <p 
corresponds to certain foliations by holomorphic discs arising from a solution of a 
certain homogeneous complex Monge-Ampere (HCMA) equation. To describe this 
case, let tt2 : D x M — » M denote the projection onto the second factor and consider 
the HCMA equation, 

{tt^uj + v^IdcVr 4 " 1 = 0, on D x M, (12) 
(vr^ + V=ld3<p)\ {t}xM > 0, ViGA (13) 
ip = i/i, on 3D x M. (14) 
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One may show that this HCMA is the Euler-Lagrange equation of an infinite- 
dimensional version of a Wess-Zumino-Witten model, given by the energy functional 



E? zw (b) = 1 -f D \vbf + j z e, 



on the space of maps b G C°°(D, G c /G), where a : dD — > G c /G is a fixed map a 
|D1] . The Lie bracket of G determines a 3- form 9 and Z is any cochain with boundary 
b(D) — bo(D) for some fixed reference map bo with the same boundary conditions 
t/j. The Euler-Lagrange equations for this functional are the WZW equations 

+ = 0, (15) 

in Euclidean coordinates q + \/—ls S D, and where d* maps sections of T*D <g> 
b*TG c /G to sections of b*TG c /G. Finally, when G and G c /G are replaced by 
Ham(M, u>) and TC, the Christoffel symbols are given by r(^,ry)| </ , = — ^^(VC, Vry), 
and the WZW equation is 

<P qq + Vss - 5|V(/? g | 2 - ^\Vip s \ 2 + {ip q , (p s }u v = 0. 

It is a perturbation of the usual harmonic map equation by a Poisson bracket term. 
Coming back to the toric situation and restricting to the space of torus-invariant 
Kahler potentials TC{T) C TC the functions (p q and (p s are commuting Hamiltonians 
and hence the WZW equation reduces to the harmonic map equation. The finite- 
dimensional WZW equation on GL(dk + 1,C) /U(dk + 1) may be written similarly 

T- l T qq + T' l T ss - (T- 1 ^) 2 - (r- 1 ^) 2 + T^Ts] = 0. 

The torus-invariance then corresponds to restriction to diagonal matrices and again 
the last term vanishes and the equation reduces to the harmonic map equation. 
Geometrically, the curvature of TC comes from the Poisson bracket and when we 
restrict to the flat subspace TC(T) the noncommutativity disappears. 

In the case of the unit disc, the normal derivative of the Green kernel is the 
Poisson kernel, whose restriction to D x dD takes the form 

Pfre^ 19 , e^ 7 ) = P r (6 - 7) = l ~f . 

v ; n 27rl-2rcos(6»-7) + r 2 

(our convention is that the Green function be nonnegative, as explained in Section 
[3]). Then we have the following more explicit statement of Theorem 1.1: 

Corollary 1.2. Let (M,L,lu) be a polarized toric Kahler manifold. Let if be a 
solution of the HCMA equation nty) -(T4\) with ip : S 1 — > TC(T) a smooth map, and 
let ifk '■ N — > Bk{T) be given on the open orbit by 

(^(re^V) = -log V \Xa{z)\ 2 h k exp ( / P r {6 - 7) log ||xa||^ d$). 

k rirzL™ K JdD 1 



aekPa 
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Then lim^— >oo 

tp k = ip in the C 2 {D x M) topology. 

The proof of Theorem 1.1 builds upon the machinery developed by Song-Zelditch 
for the study of geodesies in 7i(T). In the geodesic case, i.e., ./V = [0,1], the 
approximating Bergman Kahler potentials take the form 

= Jlog £ I^W^e-^^l^ll^)-*^''^'^). (16) 

a<=kPnZ m 



We see that the straight line segment in the case N = [0, 1] is replaced by the 
harmonic extension of the boundary L 2 norming constants in the general case. Aside 
from justifying the general formula, we need to modify the estimates to apply to 
harmonic functions on N rather than linear functions on [0, 1]. Using the localization 
lemma and the asymptotics of the peak values proved in }SoZ2j (which we recall 
in Section [2]), the uniform convergence in C 2 reduces to a verification of orders 
of amplitudes where the analysis is carried out separately in the interior of the 
polytope and near its boundary. Since many details are similar to the geodesic case, 
we concentrate here only on the novel features and the reader of Section [5] would 
benefit from some familiarity with |SoZ2] , The reader is also referred to [R] . Chapter 
3, where the results of the present article appear with greater detail. 

Let us make note of one more relation between the geodesic segment problem and 
the harmonic mapping problem. In both cases a key aspect of the toric situation 
is that the Legendre transform linearizes the harmonic map equation. This was 
known previously for geodesies iGl |S2] (see also |SoZ2j for a simple proof), but 
is observed for the first time here for general harmonic maps. We refer the reader to 
Section [3] where we also observe a generalization (|36p of a well known formula from 
convex analysis and show that the Eells-Sampson harmonic map flow is Legendre 
transformed to the usual heat flow. This fact is quite a remarkable property of 
the toric situation and does not hold for general variational problems. It follows 
that one can explicitly solve the WZW and harmonic map equations in terms of 
the associated symplectic potentials. We make crucial use of this in proving the 
convergence of the Bergman harmonic maps, and it is the reason why we do not 
founder amid regularity problems as in the general projective case. Our results give 
the first proof of convergence for higher dimensional harmonic maps into spaces of 
Kahler metrics. We hope to discuss in a separate article convergence results for 
WZW maps on general projective manifolds. 

Y.A.R. would like to thank Gang Tian for suggesting this problem to him and 
for his advice and warm encouragement, and Bo'az Klartag for enjoyable lectures 
and discussions on convex analysis. Both authors would like to thank the Technion, 
where this work was begun, for its hospitality in Summer 2007. This material is 
based upon work supported in part under a National Science Foundation Graduate 
Research Fellowship and grant DMS-0603850. 
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2. Background results 

We begin by recalling some basic facts regarding toric varieties relevant to our 
setting. We refer the reader to [Cal ISTZ1|, [SoZ2] for more details. 

We will work with coordinates on the open dense orbit of the complex torus given 
by z = eP^+v^ 9 , with p,9 G M. m x (S 1 )" 1 . The real torus T 9* (S 1 )" 1 acts in a 
Hamiltonian fashion with respect to u. The image of its moment map [i is a Delzant 
polytope P C M m . Let x denote the Euclidean coordinate on P. The polytope is 
given by P = {x G M m : l r (x) := (x,v r ) — X r < 0, r = 1, ...,d} where v r is an 
outward pointing normal to the r-th (m — l)-dimensional face of P (also called a 
facet) and is a primitive element of the lattice Z m . 

The toric monomials {Xa(z) '■= z a } a& kPnz m are an orthogonal basis of H°(M, L k ) 
with respect to any element of Bk(T). Hence a toric inner product, equivalently a 
point in Bk{T), is completely determined by the L 2 norms (up to k n /V), or norming 
constants, of the toric monomials — 



Qh*(a) := llXallh* = ! \z a \ 2 e- k *dV h 
J(C*) m 



(c*y> 

Unlike in the Introduction here we let h = e _</3 with ip a local Kahler potential on 
the open orbit (that does not extend globally). Define the normalized norms of the 
monomials 



V h k(a,z) 



\Xa(z)\ 2 



ii k 



Iv II 2 
\Xa\\ h k 



and their peak values 



V hk {a) := ^J(k))\l\ (17) 

I \ Xa | \j l k 

In order to complete the proof of Theorem 1.1 we will need some of the tools de- 
veloped by Song-Zelditch, that we now recall. First, an asymptotic expression for 
Vhk (a) for families of toric Bergman metrics. This expression is sensitive to the dis- 
tance of a/k to the boundary of the polytope. Recall that the k-th Bargmann-Fock 
model on (C, ^J—ldz A dz) is given by the holomorphic functions that are 1? with 
respect to the Hermitian metric /ig F = e _fc ' 2 ' and a basis is given by all monomials 
z a with a £ Z+. One may compute that 

p k ( a ) = k e - a ?-r. (18) 

BF OL\ 

Let 5k = k - Denote by F& k (x) = {r : l r (x) < 5k} the index set for those 



facets to which x is — ^-close, and let 5 k (x) denote the cardinality of this set. 



yk log k 

Set 
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where 5 ip (x) is defined in (f26l) below and put 

These two terms are the far and near contributions to the asymptotics of the peak 
values V h k(a): 

Lemma 2.1. (See |SoZ2| . Propositions 6.1, 6.5.) Let 5 k = k ■ Let {h t }teK be 
a family of metrics with K compact. Then there exist C > independent of t such 
that for any 5 S (0, 

V h ,(a) =Ck^ m - 6 ^^^)V BF , 5k {^ + M^h t )\ (19) 

where R k = 0{k 6 ^^). This expansion is uniform in t and may be differentiated 
twice to give for j = 1,2 and for some amplitudes Sj of order zero the expansion 

(|)\*(«) = OnA^-^y^PWf )(S 3 (t,a,k) + i? fc (f M). (20) 

Next, recall the following asymptotic localization of sums result: 

Lemma 2.2. (See [SoZ2] . Lemma 1.2.) Let B k (y,a) : kP n Z m -> C be a family 
of lattice point functions satisfying \B k (y, a)\ < Cok M for some Co,M > 0. Fix 
S £ (0, 1/2). Then there exists C > such that 

B k (y,a)V hky (a,z)= B k (y, a)V h * (a, z) + 0(£T C ). 

kPnZm a:\f-n y (z)\<k S -h 

3. LEGENDRE TRANSFORM TO HARMONIC FUNCTIONS AND LEGENDRE DUALITY 

OF GEOMETRIC FLOWS 

Over the open orbit of M, a toric Kahler potential may be identified with a 
convex function on W 1 , ip = <p(p) in the logarithmic coordinates of Jj2j By abuse 
of notation we will frequently identify Ti.(T) with the local torus-invariant Kahler 
potentials defined on the open orbit. The gradient of ip(p) = tp{e p ) is a one-to-one 
map, identified with the moment map p, whose image is P. The Legendre transform 
appeared in the context of symplectic toric manifolds in the work of Guillemin |Gu| 
and this tool lies at the heart of this section (some of our work in this section is 
also related to work of Semmes [SI]). It takes Kahler potentials <p on the open 
orbit to symplectic potentials C(p = u v , that are defined as convex functions on P 
with logarithmic singularities on dP, and relates the moment map, local symplectic 
potential u^(x) = u !fi (p(p)) and local Kahler potential as follows: 

u v (x) = {x,2\o^ > p- 1 (x)) - ip(p~ 1 (x)) 

(21) 

= (x,(V^)- 1 (x))-^((V^)- 1 (^))- 
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In addition 

(Vv?) _1 (x) = Vit(x), (22) 

and 

(VV)~ 1 |(Vv)- 1 (*) =v2u l*- ^ 23 ) 
Any symplectic potential u can be written as Uo + f, with respect to the canonical 
potential 

d 

u (x) = J2 l k{x)logl k (x) (24) 

k=l 

introduced by Guillemin, with / smooth up to the boundary [Gu] . Just as for Kahler 
potentials we may define the space of global symplectic potentials: 

CHiT) = { f G C°°{P) : n + / = dp with if g H(T)}. (25) 

We will sometimes, by abuse of notation, identify elements of this space with their 
local symplectic potentials in the same manner as with Ti(T) itself. 
Letting G tp (x) = V 2 u ip (x) we have the following formula of Abreu 

detG" 1 = 6 v (x)-Ii$ = Mx), (26) 

for some positive smooth function 5^ [A] . 

Let n = diniR N and denote by y\ , . . . , y n local coordinates over some coordinate 
patch U C N. We assume that N is oriented as a manifold with boundary, i.e., 
that the orientation on dN is the one induced from N. Recall that there always 
exists a Dirichlet Green function G(y,q) G C°°(N x N \ diag(A r )) for the Laplacian 
An ■= Ajvj on such a manifold [Aul IGTj . If v € C°°(dN), the equations 

Anu = 0, on N, 
u = v , on dN, 

have a unique smooth solution 

«(■)=-/ v{l)du{ q )G{ ■ , q)dV 9NJ (q) (27) 

(our convention will be that G(y, q) is positive in the interior and vanishes when 
q is in the boundary), with u(q) = i /i (q)-^ji\q an outward unit normal to dN in 
N (coming from a Riemannian splitting TN\qn = TON © Nqn, where Nqn is the 
normal bundle to dN in N), and where we let d ll r q )G(y,q) := u(q)G(y,q) < be 
the normal derivative with respect to the second argument. 

Let T c ah denote the Christoffel symbols of (N, f) with respect to local coordinates 
y 1 , . . . ,y n . Recall the following expression for the Christoffel symbols of (TL,gi2). 
Our proof is a slight variation on those in [CCht ID1| . 

Lemma 3.1. For every e, f G T V H we have T(e,f)\ v , = — ^^(Ve, V/). 
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Proof. Recall that the proof of Koszul formula for the Levi-Civita connection of a 
finite-dimensional manifold [Pj, page 122, carries over to infinite-dimensions to show 
that if a Levi-Civita connection exists it is unique. Regard the functions c, e, / as 
constant vector fields on Ti and let D denote the Levi-Civita connection of g^i. 
Therefore the corresponding brackets vanish and we have 



Since 



we have 



2g L 2(D c e,f)\ tp = cg L 2(e,f) - f g L i(e,c) + eg L 2(f, c). 

at ov j M v j M 

2g L 2{D c eJ% = i I (efA v c-ceA v f + fcA v e)co n 

v JM 

= 1/ (efA v c- fA v (ce) + eA v (fc))u; n . 

v JM 

g v (Vc,Ve)fu n . 



(28) 



(29) 



1 

V -I M 



It follows that 



D c e = -\g v (Vc,Ve) 



Finally, this expression is symmetric hence D is torsion-free, and it is also compatible 
with g L 2 since / g L 2 (c, e) = g L 2 (D f c, e) + g(c, D f e) is just 



l - I ec\A g Ju n = -±- f I( e V/-Vc + cV/-VeK\ 

v JM v JM 



□ 



Several authors observed previously that the Legendre transform linearizes the 
geodesic equation, so that a geodesic ipt with endpoints 930 ,fi is given by tp t = 
£ _1 (£(^o + t(C(pi — Cipo)) [DH [Gl IS2] . We now observe that under the Legendre 
transform, a harmonic map into 7i(T) is mapped to a family of symplectic potentials 
that are harmonic functions in the N variables. The key point is that the Legendre 
transform eliminates the Christoffel symbols in a variational sense. 

Proposition 3.2. Let ip : dN — > 7i(T) be a smooth map. There exists a unique 
harmonic map ip from N to TL{T) that agrees with tp on dN. Moreover, ip = C~ 1 u 
where u £ C co (N x P\ dP) satisfies Ajyu = and u\qn = Cip. 



Proof. The proof of the one-dimensional case ( |SoZ2j . Proposition 2.1) carries over 
with minor changes. Indeed, harmonic maps into TC(T) are stationary points of the 
functional, 



E(ip) 



N 



\dip\ 2 alV NJ 



NxM 



jab di P dl P 

dy a dy 



(30) 
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First considering a variation of (I2ip at p = p(x) = (V92) (x) yields 



du 

QyCl 



da 
dy a 



n 

i=i 



dpj dip 
' dy a dy a 



dip dpj 
^ dpj dy a 



dip 

QyCL 



since V(p{p) = x. 

Next, torus invariance allows us to integrate instead over the polytope and we 
have using (22) and (23) that 

(Vp)*(a;£) = (Vp)*((det VV)dpi A d(9i A • • • A dp m A d0 m ) 



(Vu); 1 ((det X7 2 ip)dp 1 A d#i A • • • A dp m A cZ0 m ) 
(V«)*((det V 2 ip)dpi A d6»i A ■ • • A c2p m A <i0 m ) 
(det VV)(det V 2 ^)^ 1 A • • • A fix™ = <ix. 



(31) 



This means that the metric g^2 is pushed-forward to the 'Euclidean' metric on 
CH{T). Therefore the functional E(ip) equals 



/ 

JN 



f 



ab 



du du 



-dx A cLVnj 



\du\ 2 dViy,f, 



(32) 



'NxP dy a dy b j N 
which is equal to the energy of the map u : iV — > C7i{T). Since the target space is 
now flat with vanishing Christoffel symbols the Euler-Lagrange equation is A^u = 0. 

Finally, note that since u\qn = CiP\qn is convex on the boundary it is also convex 
in the interior of iV: Observe that from (27) it follows that the Hessian of u (in the 
P variables, namely x) for every y £ N is given by 



V 2 u{y,x) 



dN 



V 2 u(q,x)d u{q) G(y,q)dV dNJ (q), 



and since —d l/ f q \G(y,q) > (see the paragraph after (27) ) it follows that \7 2 u(y,x) 
is therefore a positive-definite matrix. Therefore ip := C~ 1 u is in TL(T) and solves 
the harmonic map equation with boundary values tp, as required. □ 

The Eells-Sampson harmonic map heat flow |ESj on the space of smooth maps 
from (JV, /) to (H(T),g L 2) is given by 

dtp = f ab d r d yb <p - f ab T c ab d y ^ - y ab g(Vd y aip, V d yh ip) , (33) 
while the heat flow on the space of symplectic potentials CTC(T) is given by 



d t u = A N u. (34) 

Note that the equations (33) and (|34p hold without change for the global Kahler 
and symplectic potentials, respectively. 

We record the following result although we will not make use of it for the proof 
of the main theorem. 
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Theorem 3.3. Under the Legendre tranform the Eells-Sampson harmonic map flow 
k°3°J\j on the space of Kahler potentials 7i(T) is mapped to the heat flow fS^l) on the 
space of symplectic potentials CH(T). 



Proof. As above, taking a variation of (pTj) yields 



du 
~dt 



dip 
"dt 



p{x) 



(35) 



Intuitively, the equality of the energy functionals (1301) and (|32j) then suggests that 
their Euler-Lagrange equations should coincide, however up to a sign, coming from 
the fact that an infinitesimal variation 5ip(p) in one corresponds to an infinitesimal 
variation —5ip(x) in the second. More precisely, we have 



- A N u = f ab d y ad ybV - f ah T c ab d y oip - y ab g(Vd yal p, Vd yb <p). (36) 

To demonstrate (I36p . recall first the following formula for the second variation of 
the Legendre duals of a family of convex functions (parametrized by t, say) that 
have the same gradient image, that follows by taking a variation of (|2ip and using 

mi- 



d 2 u 



d 2 y 



dt 2 
d 2 ip 



dt 2 
d 2 ip 



dt 2 



or more succinctly 



(V(p)-l(aO 



E 



SV d({V?)-\x))j 



J dtdpj 
l J 



dt 



(Vvp)- 1 ^) ^ dtd Pj 



E 



d 2 ip d(du/dxj) 



Dt 



(37) 



(y<p)-Hx) 



(V(dtp/dt)\ {Vv}) - Hx) ,V(du/dt)\ x ), 



Cp = u+ (Vip, Vm). 



(38) 



Now, the terms that are linear in the first derivatives on each side of (|36p are equal 
to each other by the first variation formula for the Legendre transform (|35p . Next, 
fix y G N and choose coordinates on N for which f ab = 5 ab at y € N. Then (|38p 
gives 



-dyCidyaU = dyadyblf + (VdydU, V 8yO. If) . 



But now 



d x .{d y *u{x)) = -d Xj (d y aip(p)) = -d Pk d y aip • d Xj p k = -d Pk d y aip • d Xj (Vu) k . 

Therefore using (|23|) and the fact that is represented in coordinates on the open 
orbit by V 2 (/? we see that (|36p holds. Thus, the Legendre transform sends solutions 
of d33j) to solutions of (JMD- □ 



BERGMAN APPROXIMATIONS OF HARMONIC MAPS 



13 



In general one does not expect the Euler-Lagrange equations of two equal func- 
tionals defined on two different spaces to transform to each other. In our situation 
this does happen and in essence is due to the fact that the Legendre transform 
eliminates the Christoffel symbols not only in a variational sense but pointwise. 

Observe that Equation ([36]) generalizes the well-known formula ([38]) from convex 
analysis for the second variation of a family of convex functions on M m parametrized 
by (R, dx) that have the same gradient image. The factor ^ in our formula comes 
from the conventions we used to relate the Riemannian and Kahler metrics. 

4. The approximating sequence 

Given a harmonic map ip : N — > 7i{T) we now define the purported approximating 
sequence of harmonic maps : N — > Bk(T). First, given a family of toric Kahler 
metrics ip parametrized by dN we 'project' the family pointwise by FSfc oHilbfc onto 
Bk{T) to obtain a family of toric Bergman metrics parameterized by dN. Each 
of these metrics is determined by its I? norming constants, hence by the diagonal 
matrices 

diag(Q h k (a)) aefcPnZ m, q £ dN. 

4>(q) 

For each a, we solve the boundary problem 

AA a (y) = 0, y G N, 
K = 1ogQ hk (a), qedN. 

We then map back to via FSk to obtain the family 

MV,*) = l l °S E e ' Xaiy) \Xa(z)\ 2 hk G B k {T) 

a£kPDZ m 

of harmonic maps alluded to in Theorem 1.1. This may be written somewhat more 
explicitly in terms of the Green kernel: 

<Pk(y,*) = riog \Xa(z)\ 2 h k exp ( J d u{q) G{y,q)\og\\xa\\l^ q dV dNJ {q)^. 

Our first aim is to prove the C° convergence by showing 
<Pk(y,z) - <p(y,z) = 

ilog \Xa(z)\ 2 hk exp(/ d Hg) G(y,q)log\\xa\\l,,JV dN j(q)) =0( 1 ^). 

We begin by rewriting the sum in a convenient way. 
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r Q h k (a) 

K k (y, a) := exp ( - / d u r q) G(y, q) log v(y) , dV dNJ (q) 

V JdN ^h k , A a ) 



Put 

,(«) 

(40) 

hi / 
Then proving (139p is equivalent to proving 

-log T n k {y,a)V h u i (a,z) = 0{\ogk/k). (41) 

aekPnz™ 

Put Uy := u v fy\ = u(y, •), for y £ N. In light of the results in the geodesic case 
[SoZ2j we expect the asymptote of IZ^ to be the following: 

Definition: Let the metric volume ratio be the function on N x P defined by 

(f det ( x^] \ 

- \ J m d Hq) G(y, Q ) log detV , u y q{x) dV d Nj( q )) • 

Note that by (|26|) we have 

i i f dj,i q \(x) 
Koc(y,x) = exp ( - 2 / d u{q) G(y,q) log ^dVajv./fe)), 

and therefore IZoc G C°°(N x P) (up to the boundary). 

In light of Lemma 12.11 it will be useful to express the ratio IZk in terms of the 
functions V h k (a) (117j) in the following form: 

Lemma 4.1. One has 



^-fc(y, a) = exp ( - / d u{q) G(y, q) log - — q ——dV dNJ {q) 
Proof. By definition, 

/ /■ Qh k ( Q ) \ 

7£ & (i, a) = exp I / d v i q) G(y, q) log ^ w . rfVajv,/(g) ) . 

Specializing (|21|) to the lattice point a we have 

u<p(a) = (a, 2 log iT 1 (a)) - ip(fj,' 1 (a)), 

implying that 

log Q h k {a)V h k (a) = ku(a/k). (42) 
Since u is harmonic in y it follows that 

log Q h k (a)V h k (a) = - / d v{q) G{y, q) log (a)V h k (a)dV dN j(q), 
which together with the definition concludes the proof. □ 
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5. Proof of Theorem 1.1 

Recall that the usual elliptic regularity theory applies to the operator A^r [A-n\ 
IGTj . Namely, there exists C = C(N,f) such that the Schauder estimates hold: 

\\ u \\c^/2(N) < C(\\ u \\c°(N) + \\v\\c^/2(dN))- ( 43 ) 

Moreover, the maximum principle implies that ||ii||c°(.2V) < IMIc (aw)> an d so the 
estimates are only in terms of ||i>||c*2,i/2(gjvy 

First, we need the following asymptotic regularity for the coefficients 1Zk(t,a). 

Lemma 5.1. There exists a positive constant C > such that for all k, y, a one has 

1/C <TZ k (y,a) <C. (44) 

Moreover, fixing 5 E (0, 1/2) one has 

n k {y,a)-lZ 00 {y^/k)=0{k & - 1 2) (45) 

uniformly in (y,a), and log1Zk(y,a) is uniformly bounded in C 2 (N). 

Proof. The Bargmann-Fock terms in (|19|) depend only on the geometry of P and 
not on y E N and so are cancelled in the ratio 7Zk(y, a). Therefore, by Lemma l4.ll 

\ogK k (y,a) = f d N -d u{q) G(y, q) log V h k dV dN j(q) - log V h h 

y>W <p(y) 

= \ IdN - d ^) G iy^) lo ^{^{Q)^)) dV dNj{q) ~ 5 lo g(^(y)(D) 

+ IdN ~ d u(q)G(y, q) log (l + i?fc(§ , h^ {q) ))dV dNJ {q) 
-log(l + J R fe (s hv{y) )). 

The first two terms simplify to 

- / -d v{q) G{y,q)\og-^^dV dNJ {q) = logK^a/k), 

2 JdN d *P{q)\k) 

and so the Schauder estimates (f4*3j) together with Lemma 12.11 imply Equation (|4"5j) 
and hence also the uniform estimate (|44p . 

We now turn to prove the higher derivative estimates. A first derivative of the 
fourth term yields, according to (i20|) . 

Si(y,a,k) + R k (^,h^ y )) 

l + R k (j i 

and this is uniformly bounded according to Lemma fe.ll In a similar fashion it follows 
that second derivatives are uniformly bounded as well. Finally, the Schauder esti- 
mates (|43p may be invoked again for the third term and these will be uniform since 
the same argument as for the fourth term implies that || log(l + i?fc(^, h^^))\ \c 2 (8N) 
is uniformly bounded. □ 
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Note that the estimate (I44D immediately implies the C° convergence of (fk to ip 
with remainder as in (|39|) since we have a asymptotic expansion for the Szego kernel 
that to first order equals 

U hk (z,z)= Y, r hk (a,z) = l + 0(k- 1 ). (46) 

v\v) * — ' <p(y) 
Q6fcPnZ m 

We now turn to showing C 1 and C 2 convergence. In other words, our aim is now 
to show that the C 2 (N x M) norm of the left hand side of flU} is still 0{k 5 ~^). In 
order to prove these estimates it is crucial to make use of some cancellations. These 
can be understood as follows. When one replaces all the coefficients 1Zk(y,a) by 
a constant, one reduces to the case of a zero-dimensional map, or equivalently to 
the known asymptotic expansion of the Szego kernel that may be differentiated any 
number of times with a small error. Now there are two cases. When a coefficient 

T^kiy, ot) or a derivative thereof only multiplies a normalized monomial V h k (a, z) 

•p(v) 

it is enough to use the uniform estimates given by Lemma I5.ll and one does not 
need to keep track of error terms. However, as is usually the case, if the coefficient 
lZk(y,a) or a derivative thereof multiplies another term that itself depends on k, 
one needs to keep track of the remainder of order 0(k s ~2) given by Lemma I5.ll 
When such an error is introduced we simultaneously apply Lemma I2.2I to localize 
to those lattice points satisfying \a\ < k? +s . Remembering the overall factor of r 
one then estimates the remainders thus introduced. 

Let us now consider derivatives solely in the M-directions. A derivative of (|41 j) 
in the pj directions amounts to multiplying each coefficient in the sum (|4ip by a 
factor of 

k((V<Py)(z) ~ f )j = k(p y (z) - f ) j . 

(recall that the moment map p y is the gradient of the open orbit Kahler potential 
Namely in the interior of P one has, 

3 1 EaefcPnz™ k M z ) ~ f ot)V h u (a) 

j^-(fk - <p)(y, z) = ^ — — — — • 

d Pj k 22aekPnz^ n k{y,a)V h k (a,z) 

v(y) 

The factor of k is cancelled by the overall factor of \ and the coefficients TZk(y,C() 
are uniformly bounded due to (|44|) . Thus by Lemma [2 .21 one may restrict to those a 
such that \p y {z) — j\ < k s ~z (introducing an error 0(k~ M ) for some large M > 0). 
It follows then that 

-£-{<Pk- <p)(y,z) = 0(k 5 - 1 *). 
d Pj 

Near the boundary of P one performs the same computation but with respect to 
the slice-orbit coordinates (the same remark applies to all the computations in this 
Section). Note that the argument reduced to the one in [SoZ2j, §7.2, once we had 

(03D- 
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Next, we consider second derivatives in the M-directions. Symmetrizing sums 
(sec [SoZ2j, §8) one obtains in the interior of P, 

d 2 , \ / \ d 2 ip(y,z) 



dpidpj ' Opidpj 

1 3 £ a ,/36fcPnz™(«i ~ &)( a i - Pj)K k (y,a)ll k (y,f3)V hk (a,z)V hk ((3,z) 

_|__ v{y) w) 

k J2aekPn^ n k(y,a)V h k (a,z) 

<p(y) 

Equation (|4"5j) allows to reduce the computations to those in the case N = [0,1]: 
After localizing (keeping only those a that are 0(fc2 +<5 )-close to kfiy(z)) we have 
the estimate \ («j — fy) (pij — Pj) = 0(k 2S ). We then replace the coefficients TZ k (y,a) 
by the uniform constant (independent of a) TZ^y , p, y {z)) at the price of an error 

0(£; 35 ~2). But now what is left is then precisely cancelled by — ^0§j^ (up to 

an error of 0(k~ 2 )) due to the complete asymptotics of the Szego kernel of a single 
metric. To prove this last claim, we consider the situation of a family of Szego kernels 
parametrized by a compact manifold N, corresponding to the family of Hermitian 
metrics h y , y G N. In the toric situation this may be written explicitly as 

«£fcPnZ m n y K ' a£kPnZ m 

Then, ip(y,z) + ^ log Il^fc (z, z) = 0(k^ 1 ) (note the minus sign convention u> = 
-dd(f) has a complete asymptotic expansion, and a first space derivative gives 



U k (z,z):=U hky (z,z)= £ -q-J^= E ~q-^ 



2 



9<p(y, z ) , „,l-2\ dt p(y> z ) i i ^iogn fc (z,z) 1 x ^ e <Q ' p) 

— +0{k ) = — +- = T {U k {z,z)) > atj — -, r . 

(47) 



Similarly a second space derivative takes the form 
d 2 v(y,z) , 2 . _ 1 . , A . 2 



+ 0(AT / ) =-(n fc (z,z))- 



dpidpj k 



e < a .p> e (0,p> 



E WTrr?^ E 



*efcPnz™ ^* ^ /3<=kPn%™ 



(n^)r!i V ( RV a , e (a+ftp > 



(48) 

(by symmetrizing sums). In conclusion we have proved the claim, and hence the 
convergence of the second space derivatives. 

Therefore it remains to consider derivatives that also involve the iV-directions. 

We consider first one derivative in the iV-directions. One has 
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Since the asymptotic expansion (I46p can be differentiated and is uniform over com- 
pact families O [Z] we have 

v- e {a ' p) a i 



°(k ^) = TT — loglLfc = --^- + - — 7 : — 

v ' kdy a K(y) x ' dy a k ^ e <a ' p) 



(50) 



First note that the term 



A; v- e (a ' p> 



ekPnz m 



is of order 0(fe by Lemma 15. 11 Thus, we are left with the task of comparing the 
last term of (|50p with 

1 ZaekPnzm n k (y, a)-^£^d y a log Qhk ' y){a) 



We now localize the sums about the image of the moment map using Lemma 
introducing negligible errors (of arbitarily high order 0(k~ M )). Then we use Lemma 
15.11 to replace each occurrence of 1Zk(y,a) by TZooiv-, Hy{ z )) P ms an error of order 
0{k & ~ 1 / 2 ) (the error comes both from (|45p and the fact that we Taylor expand 
7£oo(y, n y (z)) about Hooiv, a/k) and use the fact that since we localized the sum we 
have \n y {z) — a/k\ = 0{k 5 ~ 1 / 2 )). In the terms involving 1Zoo(y, n y (z)) the factors of 
7£oo(y, n y (z)) actually cancel and so cancel with the last term of ([50]) after localizing 
the latter. 

It remains to show that the error term overall contributes 0(k s ~ 1 ^ 2 ) to the sum 
(|49p . To that end, because of the factor of 1/k it is enough to show that there exists 
a uniform constant C > independent of k such that 

\d y «\ogQ hk{y) {a)\<Ck. (51) 
Recall that the duality (I42jl of Lemma 14.11 implies 

d y a logQ ftfe (a) = kdya Uy (a/k) - d r logV h k (a). (52) 
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The second term of the right hand side is uniformly bounded by applying (I20p . To 
evaluate the first term recall that 



u y = - Uqd u ( g )G{y,q)dV dN j{q), y G N. (53) 
JON 

In terms of canonical symplectic potential uq of (I24p one may write u y = uq + f y for 
some globally smooth function f y G CH(T) on P and thus we have 

dyau y = dyaf y = -d y a I fqd u ( q \G(y, q)dVdNj{q)i V e N. (54) 

This is uniformly bounded according to the Schauder estimates. Combining the 
above the estimate ([5T]) follows. 
In sum we have shown that 

which concludes the case of a single iV-derivative. 

We now consider the case of mixed second derivatives. We will always assume 
a, (3 G kP n 1i n and so omit that from the summation notation in what follows. To 
simplify the notation further we will fix a point (y, z) G N x M and use the following 
abbreviations: 

d d 2 
d a ■= d y a = — , d ab := dyad yb = dya£)yb , 

g<«:P> 

R-a ■= fl k (y,a), Q a := Q h k (a), V a := V h k (a), V a : 



My) f(y) Qh k ( y )( a ) 

Symmetrizing sums again, it follows that 

-{fk - <p){y,z) - 



dy a dpi dy a dpi 



k 



+ t ; • (55) 



Localizing sums exchanges the term (ay — j3j) for 0{k^ +s ) up to an error of 0{k M ) 
for some large M > 0. By applying Lemma 15. II the term 

1 \ HcJm - Pi)w lo s (^) ^npvjpp 

is ©(A;" 5-1 / 2 ). Now we replace the coefficients TZk in (|55p for the lattice points a 
that remain (near n y (z)) by the uniform constant TZ OCl (y,Hy(z)) plus an error of 
order 0(k s ~ 1 ^ 2 ). Any term that does not multiply log 5^- is of order 0(/c <5 ~ 1//2 ), or 
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smaller, by using Lemma 15.11 Now there are two kinds of terms left to estimate. 
The first involve log with all TZk replaced by lZoo(y, fi y (z)). Then the coefficients 
IZood), fJ> y (z)) cancel out and we are left with the Szego kernel approximation of 

^dy^dfj ( U P *° ^(^~ 2 )) an d this cancels with — -g^rfer appearing in ([55]) . The second 

type of contributions comes from error terms of order 0(k s ~ 1 ^ 2 ) multiplying log 

By using ([52]) we may express log in terms of the global symplectic potentials 
and using in addition the fact that a and (3 are localized to a neighborhood of size 
comparable to k l l 2+& about k/j, y (z) we obtain 

d a log | <C + k\d a fy(l3/k)-d a f y {a/k)\ < C + kd\(a - (3)/k\ = 0(k^+ 5 ), 

(56) 

where C\ is the Lipschitz constant of the smooth function d a f y (as a function of N). 
By the maximum principle C\ is uniformly bounded in terms of the boundary data 
(i.e., FSfc o Hilb,t(V')) since d a f y is in fact harmonic in the iV-variables. This implies 
that the second type of contributions are of order 0(k 2S ~ 1 ), due to the overall factor 
ofl/k. 

Note that the symmetrization of the sums was crucial here. In sum, 



dy a dfP 



Finally, we consider the case of two derivatives in the ^-directions. This case is 
somewhat more involved than the previous ones and unlike in the case N = [0, 1] 
we also need to consider mixed ^-derivatives. We have 



dabiVk ~ <f)(y, z) = -d ab ip 

i ■ * (if) + * * (t I) * * {¥M 



+ 



We rewrite this as 



-2 



(57) 



d a b{Vk ~ <f)(y, z) = -d ab if + A + B + C, 



where 
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1 T, a ,0 K-a'PaKis'Pp ■ d ab log (K a ) + \d a log ( \ d b log 



Tic 



(59) 



, E Q KaVallpPp • \d a log ^<9fe log ^ 



(60) 



By differentiating (|46p we obtain similarly (analogously to the computations leading 
to (BHD), 



0(k 



-2\ 



Q/3 



(61) 

We now localize the sums to a ball of radius k? +s about k[i y (z) introducing a neg- 
ligible remainder /error and replace the two occurrences of 7£ 7 outside of the square 
brackets in (|57p as well as in the denominator by the lattice-point-independent con- 
stant TZoo(y, Hy(z)) (in what follows we refer to this operation as "replacement") 
introducing an error of order 0{k & ^^) for each replacement. First, observe that the 
replacement in the denominator is negligible. What remains to be checked is that 
the overall error introduced by replacements elsewhere is of order 0(k 5 ~^). 

In the term A these replacements introduce a term that is cancelled by substitut- 
ing the expression for —d a bf given by (|6ip into (|57p . In addition we introduce error 
terms. Let e 7 := TZj — TZoo(n y (z)) = 0(k s ~2 ) for each lattice point 7 in a localized 
sum. The highest order remainders are terms of the form ie 7 • (B\ +B2 + -B3) where 

Bl , 1 + Bl , 2 :=a ot i„ g (J-) + ia ol o g (|) ai o g (|2), 



B-2 



d ah \ogK a + i<9 a log ( ^ 
. Hp 



db log 



"Jog! 7^)<U>s 



7^ 



The errors introduced in the replacements in the term A are of the form ^e^Bi. The 
errors introduced in the replacements of the terms B and C are \e^B2 and ^e 7 .B3, 
respectively. 



First, t^B 2 



^c 7 ^2 — ~a) by Lemma l5.ll To bound B\ and B% we will use the 
estimate f|56[) for the derivatives of the norming constants Q a . First, it gives di- 
' 7> — 0(fe 25_1 ). Second, by squaring ([56]) we also obtain ie-y-Bi^ = 
0{k 35 ~2). Finally, the maximum principle also gives, similarly to the argument 



rectly that ^£7-63 
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proving (foTI) . that 



d ab log Q a < C + k\d ab f y (a/k)\ < C 2 k, 



and thus ^e-y-B^i = 0(k s 2 ). Altogether we have shown that all the remainders 

introduced by the replacements are of order 0(k 5 ~i) for some 5 £ (0, 1/2). 
Hence we have shown that 

dabi^k - <p)(y,z) = o(k 

and this concludes the proof of Theorem 1.1. 
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